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Abstract. Dynamical zeta functions provide a powerful method to analyze low 
dimensional dynamical systems when the underlying symbolic dynamics is under 
control. On the other hand even simple one dimensional maps can show an intricate 
structure of the grammar rules that may lead to a non smooth dependence of global 
observable on parameters changes. A paradigmatic example is the fractal diffusion 
coefficient arising in a simple piecewise linear one dimensional map of the real line. 
Using the Baladi-Ruelle generalization of the Milnor-Thurnston kneading determinant 
we provide the exact dynamical zeta function for such a map and compute the diffusion 
coefficient from its smallest zero. 
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Chaotic dynamic forces to move the attention from individual trajectory to 
evolution of densities. The study of the spectrum of the evolution operator, the Perron- 
Frobenius operator, is a well established method for the analytic treatment of hyperbolic 
dynamical systems. In particular dynamical averages can be extracted from the spectra 
of the generalized transfer operator: this approach was successfully used to compute 
physical quantities like i.e. the diffusion coefficient for extended hyperbolic dynamical 
systems or escape rate for open systems pQ. Relaxing the assumption of full hyperbolicity 
is the natural step forward: lot of efforts are still devoted to understand how to extend 
the assumptions of strong chaoticity to include the presence of regions of regular motion 
into a meaningful periodic orbit expansion. Along this line it was possible to show how 
to treat systems with marginal periodic orbits 0. In such a case marginality have a 
deep effect on the analytical structure of zeta functions and the system shows peculiar 
behavior like anomalous diffusion and power law decay of correlations [Sj. 

Even without relaxing the hyperbolicity assumption a strong limitation on the 
practical implementation of the periodic orbits program relies on the need of control 
over periodic orbits proliferation: the number of periodic orbits grows exponentially 
with the given by the topological entropy) and an analytical treatment 

requires the control over such a huge amount of information. Cycle expansions technique 
prevents the explosion of the theory organizing cycles into fundamental terms and 
curvature corrections (constructed with longer cycles and their shadowing by shorter 
orbits) that comes out to be a good perturbative expansion [1. . The main ingredient is 
then a good understanding of the cycles that build the fundamental and curvature terms 
that implies a good control over the symbolic dynamic of the system. Unfortunately the 
missing of such a control is the normal situation rather than the exception. Even for a 
very simple one dimensional map without complete branches such a control is typically 
impossible. As a paradigmatic example we will use a simple one dimensional linear map 
on the real line where it was shown [3] that the diffusion coefficient is a fractal function 
of the parameter defining it. Such a complex behaviour is the result of the intricacy 
of the symbolic dynamics for a generic parameter value. Even if transport properties 
of such maps were already derived in [S] any attempt to analytical organize the cycle 
expansions for the full parameter range has failed so far: at our knowledge until now 
it was only possible to identify a family of parameter values were symbolic dynamics 
become tractable and dynamical zeta function can be derived j^j. Here we show how 
it is possible to derive the diffusion coefficient for the full range of parameter values 
using an alternative formulation of the dynamical zeta functions presented by Baladi- 
Ruelle [Hj (their approach generalize the Milnor-Thurnston work [2j originally devoted 
to topological zeta function only). 

1. Invariant coordinates and kneading determinants 

For A > 2 let's define the map g\ : R — > R as a lift of the map : [0, 1] — ► R: 

g\(x + n) = §a(x) + n neZ (1) 



Fractal diffusion coefficient from dynamical zeta functions 
9a(x) 



1 A/ ^ 

-+A(x 

2 v T 



(2) 



Let's define f\(x) : [0, 1] — > [0, 1] to be the restriction on the unit interval of g\(x) i.e. 

fh{x)=g{x) modi (3) 
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Figure 1. Example of the map /a(cc) with A = 3.23. 

A bounce of initial condition for the map g(x) will spread out diffusively on the 
real line and the set of transport coefficients (for example the drift or the diffusion 
coefficient) can be computed with the help of a generalized dynamical zeta function pQ 
of the form: 

/ z n p /3a p \ 

r 1 (^/3) = n( 1 -nyn-) w 

W V |Apl J 

where the product runs over prime periodic orbit of the reduced map f[x) and the 
only quantities that enter the definition are the period of the orbit n p , the instability 

a p = nr=o f(f( x )) and the inte 

ger jumping number o~ p = g{x)—f{x) (that remembers 
the number of cells the trajectory jumps once unfolded on the real line). Cycle 
expansions technique pQ provides an organization of the infinite set of periodic orbits 
{p} that produces (for fully hyperbolic systems) a rapidly convergent power expansion 
of the dynamical zeta function Transport properties of the system are encoded into 
the behaviour of its smallest zero Zq(/3) around its value z (0) = 1. For example the 
drift J and the diffusion coefficient D can be computed by 
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For odd integer A the map consists of complete branches and the symbolic dynamic 
reduces to a full shift on an n- letter alphabet. In such a case the dynamical zeta 
function can be written in closed form and the diffusion coefficient easily computed. 
The situation is similar for even integer. For a specific set of the parameter A it 
was possible |Hj to control the symbolic dynamic and construct the dynamical zeta 
function while for the generic case we will expect an infinite number of pruning rules 
(changing discontinuously under parameter variation) and no simple expression for the 
cycle expansions organization of periodic orbits. On the other hand Milnor and Thurston 
[7j have shown how it is possible to control the symbolic dynamic of maps of the unit 
interval following the iteration of critical points. In particular they were able to relate 
the topological entropy to the determinant of a finite matrix (the kneading determinant) 
were the entries are formal power series (with coefficients determined by the kneading 
trajectories). Later Baladi and Ruelle |H| have generalized the result to constant weight 
(see also [9 and references therein for more general results). We will closely follow the 
notation of Baladi-Ruelle [E] to construct the dynamical zeta function for the family of 
maps /a {%) defined above. 

Fix A and rename for simplicity fix) := /a(x). Call ao < a\ < ... < a at the ordered 
sequence of end points of each branch (see Fig. (£[])). Define e(x) = ±1, whether f(x) is 
increasing or decreasing and t( constant weight for x G [of-i,^]. In our specific 

case the restriction of these functions on the i-interval will take the constant values: 

€i = 1 (6) 

z e^ 

U = — (7) 

where <7j is the jumping number associated to the i-th branch. 
Let's associate to each point x the address vector in Z 7 ^ 1 

a(x) = [sgn(x — Oi), sgn(x — ajv-i)] (8) 

and define the invariant coordinate of x by the formal series: 



0(x) = J2 



n—l 



i «(/») (9) 

n=0 Lfc=0 

where the product is intended equal to one if n = (the invariant coordinate 9 is single 
valued once we put e(aj) = 0). 

Defining (^(a^) = lim z ^ a ± <p(x) we compute the discontinuity vector at the critical 
points ai for i = 1, .., N — 1: 

K l (z,f3) = ^[e(al)-e(ar)] (10) 

The kneading matrix K(z, (3) is defined as the (N — 1) x (N — 1) matrix with 
Ki] i — 1, .., N — 1 as rows. Let's call A(z,/3) = detK(z,/3) the kneading determinant. 
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2. Relation between kneading determinant and dynamical zeta function 

Without entering into many details we restrict to the case under study and present the 
formula that relates the kneading determinant and the dynamical zeta function for the 
map f(x) defined above. Let' s call {p} the set of prime periodic orbits that have a 
critical point as periodic point (i.e. f nf, (ai) = ai). We note explicitly that there is at 
most a finite number of such prime orbits (for finite number of dj). It is possible to show 
that the dynamical zeta function is equal to the kneading determinant up to a rational 
function (see [S] for more general results and a proof of relation ([lip): 



A(z,P) = R(z,P) C\z,P) 



R(z, 13) = 1 - |(ei*i + e N t N ) 



m 

{p} 



t p (z,j3)Y 



(12) 



In their original work Baladi and Ruelle inserted the product over the set {p} into 
a definition of a reduced zeta function where they do not count such (repelling) cycles. 
Here we prefer to remain with the usual zeta function and factorize out explicitly the 
critical cycles. In particular formula (jllj) shows that the smallest zero of the kneading 
determinant coincides with the smallest zero of the dynamical zeta function for (3 — > 0. 



3. Fractal diffusion coefficient 



Let's start computing explicitly the kneading determinant for the map f\. It's easy to 
see that, for i — 1, .., N — 1: 

9{at) = (..., +, +, -, ....) + t i+1 e(a+) (13) 
6{af) = (...., +, -, -, ....) + ti e(a~) (14) 

where the change of sign in the component of the first addend appear at the i-th position. 
From the choice (0) we've t i+ i = eHi and so 



Ki{z,P) = (0,. .., 0,1,0,. ..,0) + ^ei 



e 2 a (a 



(15) 
(16) 

(17) 

L ~ u J 

= l-TrW (18) 

where in the last passage we've used the fact that for matrix W of the form ([To]) we've 
Tr (W n ) = (Tr W) n . Finally the kneading determinant simplify to: 

N-l 



Simplifying the notation we get that the kneading matrix is of the form: 

Kij = 1 - Wij Wy = a ib j 
and the kneading determinant can be easily computed to be 

Tr (W n 



A = det(if) = exp 



E 



(19) 



i=i 
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Due to the symmetry under x — > —x the map g(x) presents no net drift. In 
order to show it we compute the first derivative of the smallest zero of the kneading 
determinant (|5ji . Let's relabel the intervals by the associated jumping numbers: call it 
bj j = —M, 0, .., M and call a(j) the address (JBJ) of a point in the interval bj. Let's 
define the function Sj(z, (3) by the relation: 

M 

e?" 2 9(a+)=: £ Sj(z,(3)a(j) (20) 

j=-M 

In the symmetric case iterations of the critical point ao can be mapped to iterations of 
the other 'external' critical point a at and so 

M 



-(3/2 



9 (a 



(21) 



j=-M 

With the definition above the function Sj(z,/3) receives contributions only when the 
trajectory of the critical point ao traverses the interval bj. 

Each of the a(j) has 2M components (see |HJ) that, with abuse of notation, we 
can label by k = [(— M + ±), (-2 + |), |, (2 - i), .., (M - |)]. The kneading 
determinant can be rewritten as: 



A(z,(3) = 1 + ^EE^' \*fc>PW) - *A*>-P)&{-3)\ h 



(22) 



and the first derivative is then 

dA(z,(3) 



d(3 



k ) 



(23) 



2 =l,/3=0 



z=l,0=O 

= (24) 

were we've used: 

x(j) = [a(j) - a(-j)}- y(j) = [a(j) + a(-j)] (25) 

WU))k = [*0')U = — [270")]— fc (26) 

For (3 = the leading zero of the dynamical zeta function must be z(0) = 1 and 
then the kneading determinant must satisfy: 

A(l,0) = l + ^$>,(l,0)[£(j)]* 



(27) 



M 



(2? 



-M 



where in the last passage we've use the fact that J2k [%U)]k = 4j- In order to clarify 
Eq. ()28|) let's check it for the simplest case of odd integer A. In such a case the map 
consists of complete branches and ao = become a fixed point. The function Sj(z,(3) 
vanishes for all indexes but j = —M where 



z e 



-/3M" 



n=0 



A 



(29) 
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Using the fact that M — (A — l)/2 the relation (|25|) becomes 





a(i,o) = i-4(a-i)t^— 

v ' > a v ; A-1 



(30) 



Similar relations are found for even integer A were again the address of oq can be 
written in closed form due to the fact that / n (0) = 1/2 Vn > 0. 

In the same spirit it is possible to compute the diffusion coefficient for a generic 
value of A, rewriting Eq. (JHJ) as 

1 fd 2 A{z,p) /dA^PY 
~2 



D 



d(3 2 



-I- 



dz 



(31) 



2=1,(3=0 



where 



d 2 A(z,P) 



d(3 2 



dA(z,/3) 



2A 

1,0 k,j 



dsj(z,f3) 



y{j)k + 



dz 



1.0 



1 \ - 



dsj{z,p) 
dz 



dp 



1,0 
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■ 2 %{j)k 



1,0 



(32) 
(33) 



The evaluation of expression ()31|) for a parameter choice A e [2, 3] is shown in Fig. (J2J) 
and can be compared, for example, with Fig. 10 in [I] . 



Q 0.2 




Figure 2. Fractal diffusion coefficient as a function of the slope of the map f\ (x) 
computed from the smallest zero of the dynamical zeta function. The insert is a blow 
up of a part of the main figure 
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4. Conclusions 

The fractal diffusion coefficient that appears in some simple piecewise linear map is the 
consequence of the sensitivity of the topological entropy to parameter variation: even 
arbitrary small parameter change creates and destroy infinitely many periodic orbits 
and the cycle expansions approach to dynamical zeta function becomes hard to follow. 
Avoiding periodic orbits and using the Baladi-Ruelle generalization of the Milnor- 
Thurston kneading determinant we were able to construct the dynamical zeta function 
for arbitrary parameter value and compute from it the fractal diffusion coefficient. 
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